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Modeling  of  band  offsets  in  GaN  based  heterostructures 

H.  Until 

Istanbul  Technical  University,  Department  of  Physics 
Maslak  80626  Istanbul,  Turkey 

Abstract.  We  present  a tight  binding  modeling  of  the  band  offsets  of  GaN  based  heterostructures. 
The  model  considers  the  nonorthogonality  of  the  sp3  set  of  orbitals  of  adjacent  atoms  and  spin-orbit 
coupling  and  uses  the  Hartre-Fock  atomic  energies  and  interacting  matrix  element,  obtained  by 
fitting  existing  band  structures  for  bulk  materials,  to  determine  the  valence  band  energies  at  0 K and 
1 Bar,  which  are  screened  by  the  optical  dielectric  constants  of  bulk  materials  at  any  temperature, 
pressure,  strain,  and  composition.  The  model  compares  very  well  with  experiment  for  the  valence 
and  conduction  band  offsets  of  GaN  based  and  many  other  heterostructures. 


Introducdon 

The  band  offsets  appearing  at  heterointerfaces  dominate  various  device  properties  such  as 
injection  efficiency  in  heterojunction  bipolar  transistors  (HBT)  and  the  carrier  confinement 
in  modulation  doped  FETs  (MODFETs)  [1].  The  need  for  precise  and  reliable  prediction 
of  band  offsets  has  received  the  considerable  attention  of  the  solid  state  device  scientists 
and  engineers  over  the  years  (see  Ref.  [2]  for  a detailed  discussion).  In  this  article,  the 
band  offsets  at  GaN  based  heterointerfaces  are  determined  using  the  extended  tight  binding 
model  [3].  The  presented  model  includes  the  overlapping  of  hybrids  at  neighboring  bonds 
and  anti-bonds  such  that  it  cannot  simply  be  absorbed  into  a re-scaling  of  other  parameters 
to  find  the  valence  band  energies  which  are  screened  by  the  optical  dielectric  constants 
of  constituents  which  are  temperature,  pressure,  interface  strain,  and  alloy  composition 
dependent. 

1.  Extended  tight  binding  modeling  of  band  offsets 

Consider  the  formation  of  bonding  and  antibonding  states  for  a binary  semiconductor  from 
isolated  c-cation  and  a-anion  atoms.  One  uses  the  linear  combination  of  bond  orbitals  sp3 
hybrids  of  the  s-and  p-states  to  write  a bond  orbital  wave  function  as  |0hac  >=  uc|0ha  > 
+uc|0hc  > [3].  Here  |0h  >—  j(|s  > +\/3|p  >)  is  the  wavefunction  of  the  sp3  hybrid 
which  has  the  expectation  value  of  eh  =<  d>h  |H|0h  >=  (es  + 3ep)/4  and  uc  and  ua  are  the 
coefficients.  es  =<  <j>s\H\<j>s  > and  ep  —<  4>p\H\(f>p  >=  eop  + A/3  are  Hartree-Fock 
free  atomic  term  values,  where  A/3  is  the  crystal  field  splitting  energy  of  p-states  [4]  and 
eop  is  the  value  of  ep  without  the  splitting.  By  using  the  variational  principle  the  bond 
and  antibond  orbital  energies  can  be  obtained  and  minimizing  the  expectation  value  E with 
respect  to  coefficients  uc  and  ua  (3E/9uc  = 0 and  3E/3ua  = 0)  one  obtains  a set  of  two 
linear  equations  for  uc  and  ua.  The  solutions  of  this  set  of  equations  are  the  energies  of 
antibonding  and  bonding  hybrids,  obtained  from  : j=  { (Hjj  — ESjj)  = 0.  Here  H]  i = 6hc 

and  H22  = €ha  are  the  hybrid  energies,  H12  = H^  — Vsp3  is  the  interaction  matrix  and 
S12  = S*j  = Ssp3  andSn  = S22  = 1 are  the  overlap  integrals  between  orbitals  on  adjacent 
anion  and  cation  and  on  the  same  cation  or  anion.  The  interaction  between  the  hybrids  on 
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the  same  anion  and  cation,  or  between  the  bond  and  neighboring  antibonds  will  broaden 
the  hybrid  energies  into  valence  and  conduction  bands  at  k=0,  and  for  the  valence  band 
energy  is  written  as  [3]: 


Ev(r8v)  = Eb  — Emv  = 


^hc  “b  eha 


VSP3  ^sp3 


AEh 


(1) 


AEh  = 


(fhc  ~~  ZhaY  ~ 4 (€ha  + €ha)  Ssp3  Vsp3  + 4ehc£haS2p3  + 4V2p3 


1 1/2.. 


(2) 


where  AEh  is  the  hybrid  bonding  gap.  Minimizing  AEh  relative  to  lattice  constant  a 
((9AEh)/(9Vsp3)  = 0)  one  obtains  the  overlap  integral:  Ssp3  = 2Vsp3(ehA  + ehB).  where 
Vsp3  = k/a2,  with  coefficient  k 107  obtained  from  the  fitting  to  band  structure  of  bulk 
semiconductors.  The  metallization  contribution  to  the  broadening  of  bonding  energy  Eb 
into  valence  band  energy  Ev  is  [4]:  Emv  = (l/2)[(l+o,p)Vic+(l  — ap)Via]+2acV*,  where 
V]  = (ep— es)/4  characterizes  the  s-p  splitting  in  the  free  atom.  V*  = — 3(Vss+2Vsp/v/3  — 
Vpp)  is  the  additional  contribution  to  the  matrix  elements  between  nearest  neighbor  bonds. 
Since  the  uc/uc  ratio  defines  ionic  character  of  hybrid,  the  polarity  ap  and  covalencey  ac  of 
hybrid  bond  are  defined  as  [3]:  ap  = uc/(u2  + u2 ) 1 /2  and  ac  — (1  — a2)1/2.  Demanding 
that  </>hac  > is  normalized,  < <£hacl0hca  >=  u2  + u2  + 2ucucSsp3  = 1,  coefficients  uc 
and  uc  can  be  obtained  by  solving  the  secular  equations  using  the  bonding  energy  E = Eb 
giving  uc  and  uc  for  the  bonding  hybrid: 


(VSp3  - Ssp3  E) 


sp-3  ^sp- 

USp3 


(E  - €hc) 


Usp3 


USP3 


(VSP3  “ Ssp3E)2  + (ehc  - E)2  - 2Ssp3  (E  - €hc)  (Vsp3  - Ssp3E) 


,1/2 


(3) 

(4) 


The  band  offset  at  A/B  heterointerface  is  then  found  from  the  following  equation 


AEV  = 


(5) 


where  e oo  the  optical  dielectric  constant  of  the  constituent  semiconductors.  The  conduction 
band  offsets  at  high  symmetry  points  I V)C,  L3)C,  and  X6C  can  be  obtained  from  AEcj  = 
AEgi  — AEV,  where  AEg;  ((with  i=r,  L,  X))  is  the  bandgap  difference.  The  conduction  and 
valence  band  energies,  defined  as  the  standard  state  chemical  potentials,  at  any  temperature 
and  pressure  can  be  written  as  [5]: 


Ec(T,P) 

EV(T,P) 


Ec(0,  P0)  + C(:pT(l  — InT)  — -E 


Ev(0,  P0)  + CyPT(l  — InT)  + 


B 


P 

[p 


p2 

2B 

P2 

2B 


(1  +B')P3' 
6B2  _ 

(1  +B')P3' 
6B2 


(6) 

(7) 


where  Po=l  Bar  and  P is  the  applied  pressure.  a<;  = — B(9EC/9P)  and  av  = B(9EV/9P) 
are  deformation  potentials  of  Ec  and  Ev,  with  bulk  modulus  B and  its  derivative  B'  = 
9B/9P.  Cj?p  = C|jp  — Cpp  = — CpP  + ACp  and  C®p  = CpP  standard  heat  capacities 
of  conduction  electrons  and  valence  holes,  where  Cjp  = Cjjp  = CpP  = (5/2)ks  and 
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C’[jp  = C®p  + CpP  — ACp,  with  ACp  being  the  heat  capacity  of  reaction  which  is  obtained 
from  fitting  of  the  bandgap  Egj  (T,  P)  = Ec;  (T,  P)  — Ev  (T,  P)  to  its  corresponding  measured 
value  [6].  The  effects  of  interface  strain  on  AEC  and  AEV  can  be  found  by  substituting 
P = — 2BfCfefx  for  heterolayer  and  P = — 3Bses  = 0 for  substrate  for  a (001)  growth, 
written  as  [3]: 

AEV(T,  €)  - AEV(T)  - 2avfCf€fx  [1  + Cfefx]  - 2avfCfe^  C2fe}±  (8) 

AEci  (T,  e)  = AEcj (T)  + 2acfCfefx  [1  + Cfefx]  - 2acfCf^  2 ^ | B| ^ C2fej±  (9) 

where  efx  = (afx  — af))/af.  Here  af||  = as  and  a/±  = a f[\  — Df(af\\  — cif)/cif]  and 
Df  = 2C12/C11  with  Ci  1 and  C12  the  elastic  constants. 

The  valence  band  offset  at  AXB|_XC/BC  heterointerface  is  obtained  by  taking  anion 
energy  as  epa  (C)  and  cation  energy  as  epc(x)  = x€pc  ( A)  + ( 1 — x)epc  (B)  and  lattice  constant 
a(x)  = xa(AC)  + (l  — x)a(BC),  in  calculating  Ev(x)  of  Ax  Bi_xC  ternary.  The  composition 
effects  on  conduction  band  offset  AEc(x)  can  be  obtained  with  the  use  ofVegard’s  rule  and 
Kane’s  k • p model  for  semiconductors  [1,  7]:  (l/mn(x))  = 1 + (P2(x)/3)[(2/Eg(x))  + 
l/(Eg(x)  + A(x))],  where  mn(x)  = mnAmnB/(™nB  + (l  — x)mnA)  is  the  ternary  electron 
effective  mass  obtained  from  Vegard’s  rule  [1,  7],  mnA  and  mne  are  the  electron  effective 
masses  of  binaries  A and  B.  P(x)  = xPa  + (1  — x)Pb  is  the  momentum  matrix  and 
A(x)  = xAa  + (1  — x)  Ab  is  the  spin-orbit  energy.  Pa  and  Pb  of  binaries  A and  B are 
obtained  from  Kane’s  equation  with  measured  mn(x)  and  A(x)  for  x=0  and  1 [6].  Kane’s 
equation  can  be  rewritten  in  a quadratic  form  for  bandgap  Eg p,  whose  positive  root  gives 

[3]: 

1/2 

Egr(x)  = j/y(X)~  A<X)  + X-  (A(x)  — 3/y  (x))2  + 8A(x)/y(x)  (10) 

where  y(x)  = 3(1  — mn(x))/mn(x)P2(x).  A similar  expression  can  be  found  for  Egt 
bandgap  using  equation  (31)  of  Adachi  [7].  The  indirect  gap  Egx  is  determined  from 
Egx(x)  = Egr(x)  + (EgXA  - EgpA)x  + (EgxB  - EgrB)(l  - x).  Here  EgxA,  EgxB,  EgrA, 
and  Egr  n are  the  indirect  and  direct  bandgaps  of  binaries  A and  B. 

Results  and  discussion 

The  valence  band  offsets  at  AlN/GaN,  GaN/InN  and  AlN/InN  and  GaN/GaP  heterointer- 
faces are  obtained  from  Eq.  (5)  using  the  Hartree-Fock  free  term  values  [4],  measured 
bandgaps,  and  the  spin  orbit  splitting  energies  [4,  6].  Predictions  of  this  work  (ETBl)  are 
given  in  Table  I,  compared  with  those  of  p-p  extended  tight  binding  model  (ETBl)  [8], 
self  consistent  tight  binding  model  (SCTB)  [9],  linear  mufin  thin  orbital  model  (LMTO) 
[10],  and  the  first  principles  electronic  structure  model  (FPES)  [11]  against  the  experiment 
[12,  13].  The  model  predictions  are  compatible  with  those  of  quantum  mechanical  models, 
which  are  based  on  the  elaborate  band  structure  calculations  of  heterojunction  components. 

The  model  was  also  applied  to  Ill-Nitride  based  ternary /binary  heterostructures  to  de- 
termine the  composition  effects  on  band  gaps  and  band  offsets.  The  parameters  used  for  Al- 
GaN/GaN  hetero structure  are:  mn/mo=015  and  0.25,  mp/mo=0.40  and  0.70,  As/me  V=l  1 
and  19,  and  Egr/eV=6.20  and  3.39  for  GaN  and  AIN,  respectively  [14,  15].  Although 
there  is  no  data  available  for  AEV  at  AlxGai_xN/GaN  heterointerface,  prediction  of  Eq. 
(10)  for  the  direct  bandgap  Egp  of  AlxGai_xN,  is  found  to  be  in  good  agreement  with 
experiment  14], 
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Table  1.  Comparison  of  this  work  (ETB1)  with  those  of  p-p  extended  tight  binding  (ETB2)  model 
[8],  self  consistent  tight  binding  model  (SCTB)  [9],  linear  mufin  thin  orbital  model  (LMTO)  [10], 
and  first  principles  electronic  structure  model  (FPES)  [11],  against  the  measured  valence  band 
offsets  [12,  13].  (All  values  are  in  eV.) 


System 

ETB1 

ETB2 

SCTB 

FPES 

LMTO 

Experiment 

AlN/GaN 

0.26 

0.33 

0.17 

0.84 

0.85 

0.57  ±0.22  [12] 

AlN/InN 

1.52 

1.66 

0.13 

1.04 

1.09 

1.81  ±0.20[12] 

GaN/InN 

1.27 

1.33 

-0.04 

0.26 

0.51 

1.05  ± 0.25[12] 

GaN/GaP 

1.84 

1.74 

1.89 

1.72 

— 

2.30  ±0.6  [13] 
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